NONEQUILIBRIUM FLUCTUATIONS FOR A TAGGED 
PARTICLE IN MEAN-ZERO ONE-DIMENSIONAL ZERO-RANGE 

PROCESSES 



M. D. JARA, C. LANDIM, AND S. SETHURAMAN 

Abstract. We prove a non-equilibrium functional central limit theorem for 
the position of a tagged particle in mean-zero one-dimensional zero-range pro- 
cess. The asymptotic behavior of the tagged particle is described by a sto- 
chastic differential equation governed by the solution of the hydrodynamic 
equation. 



1. Introduction 

Informally, the zero-range particle system follows a collection of dependent ran- 
dom walks on the lattice Z d where, from a vertex with k particles, one of the 
particles displaces by j with rate (g(k)/k)p(J). The function on the non-negative 
integers g : No — > M + is called the process "rate", and p(-) denotes the translation- 
invariant single particle transition probability. The name "zero-range" derives from 
the observation that, infinitcsimally, the interaction is only with respect those par- 
ticles at the particular vertex. The case when g(k) is proportional to k describes 
the situation of completely independent particles. 

The problem of the asymptotics of a distinguished, or tagged particle interacting 
with others has a long history and was even mentioned in Spitzer's paper [14] (see 
also chapters 8.1, 6. II p~5]). The main analytical difficulty is that the tagged particle 
motion is not in general Markovian due to the interaction with other particles. 
However, the intuition is that in a scale the tagged particle behaves as a random 
walk with certain "homogenized" parameters reflecting the system dynamics. 

We prove in this article a nonequilibrium invariance principle, with respect to 
a diffusion process whose coefficients depend on the hydrodynamic density, for the 
diffusively rescaled position of the tagged particle in one-dimensional zero-range 
processes when the transition probability p is finite-range and mean-zero. This 
invariance principle is the first result which captures the nonequlibrium fluctuations 
of a single, given particle in a general finite-range interacting particle system. We 
remark, however, in [4], a nonequilibrium central limit theorem was proved for a 
tagged particle in the nearest-neighbor symmetric one-dimensional simple exclusion 
model by completely different methods which rely on the special structure of the 
nearest- neighbor one-dimensional dynamics. Also, we note, in [pj, a "propagation 
of chaos" type nonequilibrium result was shown for finite-range symmetric d > 1 
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dimensional simple exclusion processes which gives the fluctuations for a tagged 
particle selected at random, or in other words the average tagged particle position; 
however, this result, which makes key use of the "averaging," does not convey the 
fluctuations of any fixed, given particle and so is weaker than the one we state in 
this paper. 

We mention also, with respect to zero-range tagged particles, previous results on 
laws of large numbers, in equilibrium [12j and non-equilibrium |10j . and equilib- 
rium central limit theorems when the jump probability p is mean-zero, ^2 jp{j) = 
[TT] , [T2] , and also when p is totally asymmetric and nearest-neighbor in d = 1 [T3] , 
and also some diffusive variance results when p has a drift YIjpU) in d = 1 
and d > 3 [13]. 

Denote by £ <E Nq , No = {0, 1, . . . }, the states of the zero-range process, so that 
£(x), x £ Z, stands for the total number of particles at site x for the configuration 



Fix an integer N > 1, scale space by N 1 and assume that the zero-range 
process rescaled diffusively, : t > 0}, starts from a local equilibrium state with 
density profile po : K — ► K+- Denote by {tt^'° : t > 0} its empirical measure. It is 
well known that w t ,0 converges in probability to the absolutely continuous measure 
p(t, u)du, where p(t, u) is the solution of a non-linear parabolic equation with initial 
condition po. 

Tag a particle initially at the origin and denote by its position at time t. 
It is relatively simple to show that the rescaled trajectory {X^ JN : < t < T} is 
tight for the uniform topology. In particular, to prove convergence, one needs only 
to characterize the limit points. 

In contrast with other models, in zero-range processes X^ is a square integrable 
martingale with a bounded quadratic variation (X N ) t given by the time integral of 
a local function of the process as seen from the tagged particle: 



where a 2 is the variance of the transition probability p(-), g(-) is the jump rate men- 
tioned before, and 77^ = t x n^ is the state of the process as seen from the tagged 
particle. Here {t x : x € Z} stands for the group of translations. In particular, if 
the rescaled position of the tagged particle x^ = X^ /N converges to some path 
x t , this process Xt inherits the martingale property from X^ . If in addition x t is 
continuous, to complete the characterization, one needs to examine the asymptotic 
behavior of its quadratic variation. 

Denote by {y„ : p > 0} the one-parameter family, indexed by the density, of 
invariant states for the process as seen from the tagged particle. Let tt^ be the 
empirical measure associated to this process: 7r^ = T X «^ t ' an d suppose that 
one can replace the local function 3(^(0))/?/^ (0) by a function of the empirical 
measure. If we assume conservation of local equilibrium for the process as seen from 
the tagged particle, this function should be h(X(s,0)), where h{p) is the expected 
value of ff(f?(0))/ry(0) under the invariant state v p and A(s, 0) is the density of 
particles around the tagged particle, i.e., the density of particles around the origin 
for the system as seen from the tagged particle. 




NONEQUILIBRIUM FLUCTUATIONS FOR A TAGGED PARTICLE 



3 



As we are assuming that /N converges to xt, since ir^ 1 = T X Nir t ' and 

tt^'° converges to p(t,u)du, we must have A(s,0) = p(s,x s ). Therefore, if the 
quadratic variation of /N converges to the quadratic variation of x t , (x)t — 
a 2 J h(p(s, x s ))ds. In particular, by the characterization of continuous martingales, 
Xt satisfies the stochastic differential equation 

dx t = ay/h(p(s,x s )) dB s , 

where p is the solution of the hydrodynamic equation, h is defined above and B is 
a Brownian motion. 

We see from this sketch that the main difficulty consists in proving the conserva- 
tion of local equilibrium around the tagged particle, without assuming any type of 
attractiveness, which is relied upon in [7j. The absence of a space average creates a 
major obstacle in this step. In contrast with the proof of the hydrodynamic limit, 
we need to replace a local function instead of a space average of translations of a 
local function. We may, therefore, only use the bonds close to the origin of the 
Dirichlet form to perform the replacement and we may not exclude large densities 
of particles close to the origin. In particular, all estimates (equivalence of ensembles 
and local central limit theorems) need to be uniform over the density. This lack of 
translation invariance confines us to one-dimension. 

The method presented here may apply to other one-dimensional mean-zero inter- 
acting particle systems. However, instead of replacing a local function by a function 
of the empirical measure, one will need to replace a current multiplied by N by a 
function of the empirical measure, as what it is done for non-gradient systems, but 
without any space average. 

2. Notation and Results 

We consider one-dimensional zero-range processes with periodic boundary con- 
ditions to avoid unnecessary technicalities. This process is a system of random 
walks on the discrete torus TV = Z/iVZ where particles interact infinitesimally 
only when they are at the same site. Fix a rate function g :No = {0, 1, . . . } — > R+ 
with g(0) = 0, g(k) > 0, k > 1, and a finite range probability measure p(-) on Z. 
The particle dynamics is described as follows. If there are k particles at a site x, 
one of these particles jumps to site y with an exponential rate (g(k)/k)p(y — x). 

For simplicity, we assume that p(-) is symmetric, but our results remain true, with 
straightforward modifications, for any irreducible, finite-range, mean-zero transition 
probability p(-). For the rate function g, we assume the next conditions: 

(LG) 3 ai > such that \g(n + 1) - g(n)\ < ai for n > , 
(M) 3 ao > 0, b > 1, such that g(n + b) - g(n) > a for n > . 

A consequence of (LG), (M) is that g is bounded between two linear slopes: There 
is a constant < a < oo such that a~ x n < g(n) < an for all n > 0. 

Denote by Qn = N N the state space and by £ the configurations of so that 
£(x), x G Tjv, stands for the number of particles in site x for the configuration £. 
The zero-range process is a continuous-time Markov chain £ t generated by 

ocjv/xo - £ 5>(*)$(e(*)) [m x ' x+z )-no] , (2.1) 
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where £ x ' v represents the configuration obtained from £ by displacing a particle 
from x to y: 

— 1 for z = x 

Z x ' y (z)= U(2/) + l for z = y 

£(z) forz^x,y. 

Now consider an initial configuration £ such that £(0) > 1. Tag one of the 
particles initially at the origin, and follow its trajectory X t jointly with the evolution 
of the process £t . Specially convenient for our purposes is to consider the process 
as seen by the tagged particle defined by r] t (x) = £ t (x + X t ). This process is again 
Markovian, now on the set Q,* N = {i] G FIn', ^(0) > 1} and generated by the operator 
L N = L e ™ + L% where L e £ v , L% are defined by 

r?(0) 



^p(z)g( V (0))A^[f(r ] ^)-m} 
77(0) 



zez 

In this formula, the translation 9 Z is defined by 

■q{x + z) for x ^ 0, — z 
{0 z r))(x) = < r)(z) + 1 fora; = 

j](0) — 1 for x = —z. 

The operator corresponds to jumps of the tagged particle, while L e ™ corre- 
sponds to jumps of the other particles, called environment. 

In order to recover the position of the tagged particle from the evolution of 
the process rj tl let N£ be the number of translations of length z up to time t: 
N* — Nf_ + 1 T] t = z r]t-. In this case, X t — J2 Z z ^t ■ As jumps arc not 

simultaneous, the processes 

Jo Vs{0) 

are orthogonal martingales and, as z p( z ) = 0, we see that X t is a martingale 
with quadratic variation 



where a 2 = ^2 Z \z\ 2 p(z). 

We now discuss the invariant measures. For each ip > 0, consider the product 
probability measures ji v = pJ^' 9 in Ojv defined by 

where g{k)\ — g(l) ■ ■ ■ g(k) for k > 1, .g(0)! = 1 and Z((p) is the normalization 
constant. Z(cp) and p v arc well defined for all ip > due to conditions (LG), (M). 
Let p = p(<p) = J r){Qi)dp, v . By conditions (LG), (M), ip p is a diffeomorphism 
from [0, oo ) into itself. Define then p p — fi v ( p ), since p corresponds to the density 
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of particles at each site. The measure {p p : p > 0} are invariant for the process £t 
(cf. |U). 

Due to the inhomogeneity introduced at the origin by the tagged particle, p, p 
is no longer invariant for the process r\t- However, a computation shows that the 
size biased measures v p defined by dv p jdp p = f)(0)/p are invariant for the process 
as seen by the tagged particle, reversible when p(-) is symmetric. Here, we take 

= ^Ooj the Dirac measure concentrated on the configuration So with exactly one 
particle at the origin, and note v p =>■ 8^ as p \ 0. 

From now on, to avoid uninteresting compactness issues, we define every process 
in a finite time interval [0, T], where T < oo is fixed. Let T be the unit torus and 
let A4 + (T) be the set of positive Radon measures in T. 

Consider the process £f —: £tN 2 > generated by N 2 Cn- Define the process 7r t ' 
mZ>([0,T],M+(T)) as 

nf'°{du) = ^(x) 6 x/N{du) , 



where 5 U is the Dirac distribution at point u 

Vo(-) 



For a continuous function po : T — > R + , define p 1 ^ , % as the product measure in 



Q.N given by p^ o ^(r](x) = k) = p pa ( x /N){ r l{ x ) = The next result is well known 
(cf. Chapter V 0; see also 0, 0). 

IV, o 

istic measure p(t,u)du, where p(t,u) is the solution of the hydrodynamic equation 

dtp = o- 2 dl<p{p) 



Theorem 2.1. For each < t < T , n t ' converges in probability to the determin- 
ation 

(2.2) 



p(0,u) = po(u). 



and <p(p) = J g{£_{0))dp 



Define now the product measure v — v N ,« in Q* N given by v N ,An{x) 



V I \ 111 ^ L i\ r t^lV Cll U V IS i 

k) = t' po ( :z ./jv)( 7 7( a; ) = k), an d let rjf^ =: n tN 2 be the process generated by N' 
and starting from the initial measure v N . Define the empirical measure ir( 
V([0,T\,M+(T))by 



(du) = — Y V^(x)5 x/N (du). 



\ I I I I I 

N 

Define also the continuous function ip : R + — > M + by 

W = /(j(#)MoK = {^ tpli 

The next theorems are the main results of this article. We first identify the 
scaling limit of the tagged particle as a diffusion process: 

Theorem 2.2. Let xf = X$ /N be the rescaled position of the tagged particle 
for the process ^ . Then, {x^ : t G [0,T]} converges in distribution in the uni- 
form topology to the diffusion {x t : t € [0, T]} defined by the stochastic differential 
equation 

dx t = a y/^(p(t,x t )) dB t , (2.3) 
where Bt is a standard Brownian motion on T. 
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Through this characterization we can describe the evolution of the empirical 
measure as seen from the tagged particle: 

Theorem 2.3. {it? : t <E [0,T]} converges in distribution on V([0,T},M+(T)) to 
the measure-valued process {p(t, u + x t )du : t € [0, T]}, where p(t, u) is the solution 
of the hydrodynamic equation \2. 2\) and Xt is given by (|2.3p . 



Recall that rj^ is distributed according to v ^ (.y Denote by F N the probability 
measure in T>([0,T],il* N ) induced by the process r)f , and by E N the expectation 
with respect to this process. Denote also by E p [h] and (h)^ the expectation of a 
function h : — * R with respect to the measure /i; when p = v p , let E p [h], (h) p 
stand for E v [h], (h)^ p . Finally, since in the next sections we consider only the 
speeded-up process r/^ we omit hereafter the superscript N. 

The plan of the paper is now the following. After some tightness estimates in 
Section [3l certain limits are established in Theorem 14. II in Section @]-with the aid 
of "global" and "local" hydrodynamics results in Sections [5] and [5]-which give the 
main Theorems 12.21 and 12.31 

3. Tightness 

To keep notation simple, in this section we assume the transition probability p(-) 
to be nearest neighbor: 

p(l) = p(-l) - 1/2 
so that a 2 = 1. Denote by C(T) the space of real continuous functions on T and by 
C 2 (T) the space of twice continuously differentiable functions on T. For a function 
G in C(T), denote by ir^ (G) the integral of G with respect to ir^: 

Trf (GO = f G[uW{du) = ^ E G(x/N)nf(x) . 

For T > 0, denote by V T = P([0, T],M + (T) x M+(T) x T x E+) the path space 
of cadlag trajectories endowed with the Skorohod topology. For N > 1, let Qm be 
the probability measure on T>t induced by the process , tt^ , x^ , (x N )t), where 
(x )t stands for the quadratic variation of the martingale . We prove in this 
section that the sequence {<5Ar:^V>l}is tight, which follows from the tightness 
of each component of (7r t ,0 ,tt^,x^, (x N ) t )- 

Let Q° N be the probability measure in P([0,T], M. + (T)) corresponding to the 
process ir^'° . As mentioned in Theorem l2.H Qf^ converges to the Dirac-^ measure 
concentrated on the path p(t,u)du, where p is the solution of (|2.2|) . Hence, the 
sequence {Q° N : > 1} is tight. 

On the other hand, as M+(T) is a metrizable space under the dual topology 
of C(T), to show that {n N : N > 1} is tight, it is enough to prove tightness of 
the projections {tt n (G) : N > 1} for a suitable set of functions G, dense in C(T). 
For G in C(T), let Qj^ be the measure in T>([0, T},M) corresponding to the process 
{vrf (G) : < t < T}. Tightness of the sequence {Q% : N > 1} follows from 
Aldous's criteria in the next lemma. 

Lemma 3.1. The sequence {Q% : N > 1} is tight if 

(i) For every t £ [0, T] and every e > 0, there exists M > such that 

su P P Ar [| 7 rf(G)| >M] <e. 
N L J 
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(ii) Let Tt be the set of stopping times bounded by T . Then, for every e > 0, 
lim limsup sup supP^ [ \ir" +e (G) - ?rf (G)\ > e 

Lemma 3.2. The sequence {Q% : N >!}, G in C 2 (T), is tight. 
Proof. An elementary computation shows that for each G in C(T) 

ft 



M t N ' G = vrf (G) - tt^G) - / ^ E {**G)(x/N)g(T, a (x))da 

■/o %(0) Jo ^ ?? s (0) 

is a martingale of quadratic variation (M N ' G ) t given by 



(3.1) 



a:eTiv\{0} 



%(o) 



[(Vjv )Z G)(0)] 2 ds 



+ jk / E^) 

+ /*EpW^7^(4 E (V^G0(^»7.W-4(VA,.fVm)i) ( /, 

•'Q »cJ ^ s *> ' »cT.. 



In these formulas, V jv>G, AnG correspond to the discrete first and second deriva- 
tives of G: 

(Vjv,,G)(u) = JV[G(« + z/JV) - G(«)] , 

(AjvG)(n) = 7V 2 EpW[ g ( u + z /^)- g ( u )] ■ 

Since the rate function g grows at most linearly and since the total number of 
particles is preserved by the dynamics, 



(M N >% = r^ p(z) ^l(l £ V N , z G{x/N)vs{x)fds + R 



N,G 
t 



where |i?f' G | < CqN 2 ^2 xe f N f]o( x ) an d Co is a finite constant which depends 
only on G, g, p and T. In particular, E N [\R^ G \] < CiN' 1 . 

Note that in contrast with the martingale associated to empirical measure 77^'°, 
due to the jumps of the tagged particle, the martingale M N ' G does not vanish in 
L 2 (P N ). In particular, we may not expect the convergence of the empirical measure 
77^ to a deterministic trajectory. 

We are now in a position to prove the lemma. Condition (i) of Lemma 13. II is a 
direct consequence of the conservation of the total number of particles. In order to 
prove condition (ii), recall the decomposition (|3.ip of tt^ (G) as an integral term 
plus a martingale. The martingale term can be estimated by Chebychev's inequality 
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and the explicit form of its quadratic variation: 



\M 



N,G 



< 



< 



N.G I 



> e 



< 



,N,G\2 



E 



N 



\n ^ 7?s(x) ) ds . 



Y {K 

2 



N,G\2 



x6Tjv 



e 2 N 



Cx\\G'\\lo6 



P0<.) 



7f £ ^ 



Cx 
e 2 N 



which converges to as N f oo and 7 4 0. The integral term can be estimated in 
the same way, using again the conservation of the total number of particles. This 
proves the lemma. □ 

It remains to consider the scaled position of the tagged particle x^ and its 
quadratic variation. We recall that x^ is a martingale with quadratic variation 



1 gMM 

V?(0) 



ds 



(3.2) 



Lemma 3.3. The process {(x N , (x )) : N > 1} is tight for the uniform topology. 
Proof. We need to show that 



lim lim sup Q 



N 



N^oo 



sup \x 

|t— s|<e 



„N 



\>s] =0 



(3.3) 



for all 8 > and a similar statement for the quadratic variation (x N )t- Recall that 
sup fc>1 g(k)/k < a < 00 and consider a symmetric random walk on the discrete 
torus T/v with jump rate a and transition probability p(-). We may couple Zf and 
Xf* in such a way that the skeleton chains are equal, i.e., that the sequence of sites 
visited by both processes are the same, and the holding times of Z N are always less 
than or equal to the holding times of X N . In particular, 



I N 

sup \x t 

It— s|<e 



..N 



< sup I z 

|t- S |<6 



.,N 



if = Zf /N. Therefore, (|3.3[) follows from the tightness in the uniform topology 
of a rescaled symmetric random walk. 

Tightness of the quadratic variation (x N )t in the uniform topology is an elemen- 
tary consequence of its explicit expression (|3.2[) and the boundedness of g (k) /k. □ 



4. Limit points and proof of Theorems 12.21 12.31 

The following, which characterizes certain limit points, is the main result of this 
section, which yields Theorems 12.21 and 12.31 

Theorem 4.1. The sequence Qn converges in the Skorohod topology to the law 
Q concentrated on trajectories { (ir® , ir t , x t , A t ) : < t < T} such that Tr°(du) = 
p(t,u)du, where p is the unique weak solution of <\2.2\ ; x t is the solution of the sto- 
chastic differential equation (|2.3[) ; Tr t (du) — p(t, xt+u)du and A t = a 1 J Q ip(p(s, x s )) 
ds. 

Proof of Theorems 12.21 and 12.31 As the limit xt is concentrated on continuous 
paths, Theorem 14. II straightforwardly implies Theorems 12.21 and 12.31 □ 
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The proof of Theorem 14.11 is now divided in a sequence of lemmatta. Denote 
by {t u : u G T} the group of translations in T acting on points, functions and 
measures. 

Lemma 4.2. All limit points Q of the sequence {Qn '■ N > 1} are concentrated 
on trajectories {(7T°, 7Tt, Xt, At) : < t < T} in which xt is a continuous square 
integrable martingale. 

Proof. Assume, without loss of generality, that Qn converges to Q. Since, by 
Lemma T3.31 {i" : N > 1} is tight for the uniform topology, Q is concentrated on 
continuous paths x t . In particular, converges in law to Xt for all < t < T. 

The martingale property is inherited by Xt because x^ converges in law to Xt 
and 

' 7M0)) 

















L Jo 



■ ds 



< aa 2 t 



»7.(0) 

uniformly in N. Therefore, Xt is a square integrable martingale relative to its 
natural filtration. □ 

Lemma 4.3. All limit points Q of the sequence {Qn '■ N > 1} are concentrated on 
trajectories {(7r°, 7r t , xt, A t ) : < t < T} in which ir°(du) = p(t,u)du, where p is 
the unique weak solution of (|2.2|) , and ir t (du) = T Xt ir®(du) = p(t,Xt + u)du. 

Proof. Assume, without loss of generality, that Qn converges to Q. The first state- 
ment follows from Theorem 12.11 On the other hand, by Lemma 14.21 and since p 
is continuous, Q is concentrated on continuous trajectories {{ir^Xt) : < t < T}. 
Hence, all finite dimensional distributions (f.d.d.) of (7r°' , x^) (and therefore 
of T x Nirf'°) converge to the f.d.d. of (n^xt) {T Xt -K^ — p(t,u + xt)du). Since 

tt^ — t x ntt^'° and since the f.d.d. characterize a measure on 2?([0,T]), the lemma 
is proved. □ 

For e > 0, denote t e = £ _1 1{(0, and a e = (2e)" 1 l{(-£, e)}(u). For I > 1 

and x £ Z, denote by r) s (x) the mean number of particles in a cube of length 21 + 1 
centered at x S Tn at time s > 0: 

\y~x\<l 

When s — 0, we drop the suffix "s" for simplicity. 

A function h : £In —> R is said to be local if it depends only on a finite number 
of sites. For a local, bounded function h : f2jv — ► K, denote by H{p) and h(p) its 
expectations with respect to v p and p p respectively. Thus, H, h : R+ — > R are the 
functions defined by 

H(p) = E Vp [h{r,)], and h(p) - [h(0] ■ (4.1) 

Also, define for / > 1 the local function Hi : Qn —> R given by 

Hi{rf) = H( v l (0)). 

Then, H l : R + -> R is the function H t (p) = E {lp [Hi], 

A local function h : jv ^ R is said to be Lipschitz if there exists a finite subset 
A of Z and a finite constant Co such that 

\h(0-h(O\ < C ^|£(x)-£'(z)| (4.2) 
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for all configurations £, £' of flN- 

Consider in particular the local function ho(rj(0)) = g(f](0))/i](0). It follows from 
assumptions (LG), (M) that ho(-) is a Lipschitz function, bounded above by a finite 
constant and below by a strictly positive constant. 

We now characterize the quadratic variation of Xt ■ 

Lemma 4.4. All limit points Q of the sequence {Qn '■ N > 1} are concentrated on 
trajectories { (ir® , 7Tf, Xt, A t ) : < t < T} such that 

r t 

2 



A t = a i>{p(s,x s ))ds 
Jo 

for allO <t <T. Moreover, A t is the quadratic variation of the martingale Xt- 

Proof. Assume, without loss of generality, that Qm converges to Q. Since (x )t is 
tight for the uniform topology by Lemma I5T51 (x N )t converges to a limit A t for all 
< t < T. By Proposition O with respect to h(r](0)) = g(r](0))/r](0) and H(p) = 
ip{p), and since for each < t < T the map it. — > f^ds J i e (x)ipi(ir s (T x a e )) dx is 
continuous for the Skorohod topology, 



lim lim lim Q 



r t 

J2 







A t - a \ ds / L t (x)tpi(n s {T x a e ))dx 



> S 



= 



for all < t < T and 5 > 0. By Lemma l4~3l ir t {du) — p(t,x t + u)du. Also, p(s, ■) 
is continuous for < s < T, and 4>i(a) — * ip{a) as I \ oo by bounded convergence. 
Then, as e | 0, e J, 0, and I | oo, we have a.s. 

L e (x) / Tpi{TT s {T x a E ))dsdx — > / ip(p(s,x s ))ds . 
Jo Jo 

It remains to show that A t corresponds to the quadratic variation of the square 

integrable martingale x t . By Corollary VI.6.6], {{xf , (x N ) t ) : < t < T} 

converges in law to {(xt,(x)t) ■ < t < T}. Since by the first part of the lemma, 

{(xf , (x N ) t ) : < t < T} converges to {(x t , A t ) : < t < T}, (x) t = A t . This 

concludes the proof of the lemma. □ 

Recall that the quadratic variation (x)t of a martingale Xt is equal to xf — x\ — 
2 J x s dx s and that (x)t can be approximated in L by the sequence of Riemannian 
sums J2j( x tj+i ~ x tj) 2 : as the mesh of a partition {tj : 1 < j < M} of the interval 
[0, t) vanishes. In particular, one can prove directly in our context the identity 
between A t and the quadratic variation (x)t- 

It follows from the characterization of continuous martingales that xt is a time- 
changed Brownian motion: 

Corollary 4.5. The rescaled position of the tagged particle {x^ : < t < T} 
converges in law to the solution of the stochastic differential equation 



dx t = a^%j}{p(t,x t )) dB t , 
where B t is a Brownian motion and p is the solution of the differential equation 



Proof of Theorem 14.11 By Section [3l the sequence Q N is tight. On the other 
hand, by Lemma l4~3l and Corollary 14. 5( the law of the the first and the third com- 
ponents of the vector (7r° , irt , Xt , A t ) are uniquely determined. Since, by Lemmatta 
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14.31 14. 4[ the distribution of the second and fourth components are characterized by 
the distribution of xt, and p(t, x), the theorem is proved. □ 

5. Global replacement lemma 

In this section, we replace the full empirical average of a local, bounded and 
Lipschitz function in terms of its density field. The proof involves only a few 
modifications of the standard hydrodynamics proof of [6, Lemma V.1.10, Lemma 
V.5.5]. 

Proposition 5.1 (Global replacement). Let r : 0,^ — > M be a local, bounded and 
Lipschitz function. Then, for every S > 0, 



lim sup lim sup '. 

E— too jV— >00 



T 1 



E 

iETjv 



N T xVsN(Vs)ds > 5 



-0, 



whe 



\v\<l 

For two measures fi, v defined on Qn (or Q* N ), denote by Ti(p\v) the entropy of 
\x with respect to v: 



TL{p\v) — sup | J fdfi — log J e^dy| , 



where the supremum is carried over all bounded continuous functions /. 



A simple computation shows that the initial entropy H{yz! t.\ \ v p) is bounded by 



CqN for some finite constant Cq depending only on po(') an d <?■ Let / t (?y) be the 
density of r/ t under P N with respect to a reference measure v„ for p > 0, and let 
fn/^\ _ ±-i r 4 fN i 



/t V ( ? ?) — * 1 Jo /i V ( 7 ?)^ s - By standard arguments (cf. Section V.2 [B]), 

H N (f t N ):=H(f t N dv p \v p )<C N and V N (f t N ) := (y/f?(-L N yffj<)) p < ^ . 

Consequently, by Chebyshev inequality, to prove Proposition 15.11 it is enough to 
show, for all finite constants C, that 

lim sup lim sup sup f V] T x V e N{'q)f{'n)dVp = 

£^0 N^oo H N (f)<CN J 

V N (f)<C/N X&lN 

where the supremum is with respect to f p -densities /. Notice that we may remove 
from the sum the integers x close to the origin, say \x\ < 2eN, because V e N is 
bounded. After removing these sites, we are essentially in the space homogeneous 
case. ProDOsition l5.ll follows from the two standard lemmatta below as in the proof 
of Lemma V.1.10]. 

Lemma 5.2 (Global 1-block estimate). 

lim sup lim sup sup f— V] T x V k (r])f(r])diy p = 0. 
k^oo n-^oo n N (f)<CN J i f~r „ 



H N (f)<CN 
T> N (f)<C/N 



\x\>2eN 
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Lemma 5.3 (Global 2-block estimate), 
lim sup lim sup lim sup sup 

fe^oo e^O N^oo H N {})<CN 
D N (f)<C/N 



hfl E jjj E \ri k {x + y)-r 1 k (z)\f{r,)dv p =Q. 



2Ne . .. 

\y\<Ne |z|>2eAr 

We now indicate the proofs of Lemmatta 15.21 and 15.31 in relation to [HI Sections 
V.4, V.5 ]. 

Proofs of Lemmatta 1 5. ,21 and \5.S\ To be brief, we discuss only the proof of Lemma 
15.21 through some modifications of the argument in [6l Section V.4] , as the proof of 
Lemma 15 .31 using the modifications for Lemma 15.21 given below, is on similar lines 
to that in [6l Section V.5]. 

In the first step of the 1-block estimate we cut-off high densities. We claim that 

lim sup lim sup lim sup sup / — T x Vk{ri)l{ri k (x) > Ajf^dvp = . 

A^oo k^oo N^ca HN(f)<CN J ^* , 

w '— \x\>2eN 

Since Vk is bounded, we may replace it by a constant and estimate the indicator by 
A~ 1 r/ k (x). After a summation by parts, the expression is easily shown to be less 
than or equal to Cq(AN)~ 1 J2 x -^ r)( x ) f° r some finite constant Cq. To conclude it 
remains to follow the proof of [SJ Lemma V.4.1], applying the entropy inequality 
with respect to v p and keeping in mind that the marginals of fj, p and v p coincide 
on sites x ^ 0. 

Define now Vk,A{v) = Vfc( J 7)l{ J T (0) — A}. By the previous argument, it is 
enough to show that for every A > 0, 

limsuplimsup sup / T7 E T ^k,A^i)f{v) dv p = ■ 

k^oo N^oo V N (f)<C/NJ iV xeJrf 

The proof is analogous to the homogeneous case. Since the origin does not 
appear, both the Dirichlet form T>n and the the measure v p coincide with the 
Dirichlet form of the space homogeneous zero-range process and the stationary 
state n p . In particular, all estimates needed involve only the functionals of the 
space- homogeneous process already considered in [5] . □ 

6. Local replacement lemma 

In this section, we replace a bounded, Lipshitz function supported at the origin 
by a function of the empirical density. 

Proposition 6.1 (Local replacement). For any bounded, Lipschitz function h : 
N — > K, and any t > 0, 



lim sup lim sup lim sup lim sup E 

l— >oo e^O e^O TV— >oo 



A' 



t 1 * 

x=l 



= 



where H{p) = E Vp [h], Hi{rf} = H( V l (0)), and Hi(p) = E Pp [Hi] . 

In the proof of this lemma, there are two difficulties. The first and the most 
important one is the absence of a spatial average, a crucial point in the standard one 
and two blocks estimates since it allows a cut-off of large densities and a reduction 
to translation-invariant densities in the estimation of the largest eigenvalue of a 
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local perturbation of the generator of the process. Without the density cut-off, 
the equivalence of ensembles, and therefore the local central limit theorem, has to 
be proved uniformly over all densities. Moreover, this absence of space average 
confines us to one-dimension. 

A second obstacle is the lack of translation invariance of the stationary state, 
turning the origin into a special site. Functions h(r)(0)) and h(ri(x)), for instance, 
have different distributions. In particular, in contrast with the original zero-range 
process, the integral J {g(i](0)) — g(rj(x))}fdi/ p cannot be estimated by the Dirichlet 
form of /. 

The proof of Proposition 16. II is divided in several steps. We start with a spectral 
gap for the evolution of the environment restricted to a finite cube. For I > 1, 
denote by A/ a cube of length 21 + 1 around the origin: A; = {—I, ■ ■ ■ , 1} and by 
L\™ the restriction of the environment part of the generator to the cube Af. 

{LffM = £ £>(y - *) 9(v(x)) [f(v x > v ) - fm 

x£Ai y£Ai 

zeAi ^ ' 

We assume above, without loss of generality, that I is larger than the range of p(-). 

Let Up 1 be the measure v p restricted to the set A;. For j > 1, denote by £a ; j the 
set of all configurations in A; with at least one particle at the origin and j particles 
in A; , and by i^a, j the measure v p l conditioned to Sa, j : 

Sa,j = {veK 1 :»7(0)>1, E^)=j}' »^A-) = ^(-l^i). (6.1) 

xeA, 

Note that j does not depend on the parameter p. 

Lemma 6.2. There exists a finite constant Co such that 

(f-J)u Al ,, < C l 2 {f(-L^f)) VAiij 

for all j > 1, all I > 1 and all functions f in ^{v^.j). In this formula, (/; f)v A } 
stands for the variance of f with respect to ^ a, j . 

Proof. This result follows from the spectral gap of the zero-range process proved in 
[8]. Since g(k)/k is bounded above and below by finite strictly positive constants, 
an elementary computation shows that 

(/;j% AiJ = inf((/-c) 2 )„ Ai ,. < a 2 if{{f -cf)^^ = a 2 (/';/% A „ 3 _ 1 

and {fi-UJ'))^^ < a 2 (f(~LT;f))^ 

provided < a^ 1 < g(k)/k < a for all k > 1. In this formula, £a ; is the generator of 
the zero-range process (|2.1[) restricted to the set A/, jUAjj-i is the canonical measure 
associated to the zero-range process restricted to the set A; with j — 1 particles, 
and f'(j]) = f(ji + Do), where Do is the configuration with exactly one particle at 
the origin and summation of configurations is performed componentwise. □ 
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6.1. Local one-block estimate. For I > 1, define the function Vi(r]) by 

vm = %(o)) - flMo)) 

where we recall h is a bounded, Lipschitz function, and H(a) — E Ua [h(rj(0))]. In 
this subsection we give the second step for the proof of Proposition 16.11 

Lemma 6.3 (One-block estimate). For every < t < T, 

ft 



N 



Km sup lim sup E 

/ — >oo N — >oo 



Proof. Since the initial entropy ^(^(.) W P ) 



Vi(rj s )ds 



is bounded by CqN, by the entropy 



inequality, 



Vi(r) s )ds 



< ^ + llogI 

7 "fN 



exp ■ 



{ 7 iv| [ V^ds^ 



where E p denotes expectation with respect to the process starting from the invariant 
measure v p . Using the elementary inequality e> x > < e x + e~ x , we can get rid of 
the absolute value in the previous integral, considering h and — h. In this case, 
by Feynman-Kac formula, the second term on the right hand side is bounded by 
(jN)~ 1 TAnj, where Ajy^ is the largest eigenvalue of N 2 L^ + jNVt- Therefore, to 
prove the lemma, it is enough to show that (jN)^ 1 Xnj vanishes, as N f oo, I f oo, 
for every 7 > 0. 

By the variational formula for \nj , 



= sup{(^/),- 7 " 1 iV(y7(-LivV7)>p} . (6-2) 



where the supremum is carried over all densities / with respect to v p . Recall that 
we denote by L e ™ the restriction of the environment part of the generator to the 
cube A ; . As the Dirichlet forms satisfy (v7 L e g° 'v7)) p < {VJ{-LnVJ))p, we 
may bound the previous expression by a similar one where Ln is replaced by L e ™ . 

Denote by /; the conditional expectation of / given {r](z) : z £ A;}. Since Vi 
depends on the configuration 77 only through {77(2) : z £ A;} and since the Dirichlet 
form is convex, the expression inside braces in (|6.2p is less than or equal to 



(6.3) 



Vihdv^ - j-'N / \Jf l {-Lf?sJf l )dv i 



where, as in (|6.ip . iA* stands for the restriction of the product measure v p to Nq 1 . 
The linear term in this formula is equal to 

j>l J 

where Vhi.j is the canonical measure defined in (|6.1|) and 

The sum starts at j = 1 because there is always a particle at the origin. Note also 
that J2j>i c l,j(f) — 1 an( i that fi.j(-) is a density with respect to v\ t j. 
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By the same reasons, the quadratic term of (|6.3p can be written as 

7" 



In view of this decomposition, (|6.2[) is bounded above by 



sup sup 

i>i / 

where the second supremum is carried over all densities with respect to v Al ,j ■ 

Recall that Vi(r)) = h(r)(0)) - H(r] l (0)). Let H t (j/2l + I) = J %(0)) dv Auj . By 
Lemma E31 below, we can replace H(rj l (0)) by Hi(rj l (Qj) in the previous expression. 
Let Vij (77) = h(r](0)) — Hi (j/2l + 1) and notice that Vij has mean zero with respect 
to v Al ,j for all j > 1. By Lemma [6.21 L A ™ 1 ' nas a spectral gap on EA,,j of order 
CqI~ 2 , uniformly in j. In particular, since h is bounded, and the inverse spectral 
gap of L A nv is order I uniformly in j, by Rayleigh expansion [6l Theorem A3. 1.1], 
for sufficiently large N, 



* l-2||Xw^ / ^D- 1 ^* 



< 



uniformly in j > 1. By the spectral gap of L e ™ again, this expression is less than 
or equal to 



because h is bounded. This proves that (16. 2| vanishes as JV | oo, / | oo, and 
therefore the lemma. □ 



Lemma 6.4. For bounded, Lipschitz function h : No 



lim sup sup 

l~>oo k>0 



E, Al , k [h(viO))]-E Vk/lAil [h(r,(0))} 







Proof. Fix e > and consider (l,k) such that fc/|A;| < e. We may subtract h(l) 
to both expectations. Since h is Lipschitz, the absolute value appearing in the 
statement of the lemma is bounded by 

C(h){ J {r,(Q)-i}du Mtk + J MO) - l}dv k/ \ H ) • (6.4) 

Note that both terms are positive because both measures are concentrated on con- 
figurations with at least one particle at the origin. We claim that each term is 
bounded by a 2 e. 

On the one hand, since v p {dr\) = {r)(0)/p}fi p (dr)), the second term inside braces 
is equal to p^ 1 J i](0){r](0) — 1} d(Xk/i Al i, where p — k/\Ki\. since k < ag(k), we may 
replace 77(0) by ag(r)(0)) and perform a change of variables rj 1 = n — do to bound 
the second term in (|6.4p by a<p(p) < a 2 p. 
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On the other hand, by the explicit formula for ^A ; .fc, the first term in (|6.4|) is 
equal to 

n^/ n^. 

xeA, av /v " x<£A t Jln " 

where both sums are performed over Sa,,^. Replacing r](0) by a ±1 g(77(0)) in the 
numerator and in the denominator, we obtain that the previous expression is less 
than or equal to a 2 E PAi k _ 1 [i](0)] < o?k/\Ki\. In last formula, /XA^fc is the product 
measure [i p conditioned on the hyperplane S A k = {£ e Nq 1 : X^eA = 
For (Z,fc) such that fc/|A;| > e, write 

^ Ai , fc [%(o))]-^ /|A(l [%(o))] = ^^.^^(o))]-^^,^^))]} , 

where p = fc/|Aj| and h'(j) = h(j)j. By Corollary 6.1 (parts a,b) [5] the last 
difference in absolute value is bounded by C(/i)e -1 / -1 . □ 

6.2. Local two-blocks estimate. In this subsection we show how to go from a 
box of size I to a box of size eN: 

Lemma 6.5 (Two-blocks estimate). Let H : M+ — » K 6e a bounded, Lipschitz 
function. For every t > 0. 

eN 



I {H(rjl(0))-±*£H( V l s (x))}ds\]=0. (6.5) 

*' 6 x=l 



lim sup lim sup lim sup E N 

l— >oo e^O iV— XX3 

The proof of this lemma is very similar to the proof of Lemma 16.31 The expec- 
tation in (|6.5p is bounded by 

' X - - C(H)(2l + l) 



^ J2 ® N [\[ {H(r, l s (0))-H(r, l s (x))}ds 



£C=22 + 1 



eA 



Following the proof of the one-block estimate, we see that it is enough to estimate, 
uniformly in 21 + 1 < x < eN, the quantity 

sup 



=U P {(Vl, x f) p - A 7 - 1 ( V ^(-iArV7)>p} I 



where the supremum, as before, is over all density functions / with J fdv p = 1 and 
Vz :2 , is defined by 

V l>x ( V ) = H(r, l (0))-H( V l (x)). 

Notice that the blocks A; and A; (x) =: {—l + x,...,l + x} are disjoint. Let L e ™ 
be the restriction of L e ^ v to the set A/ x = A; U A; (x) and define the operator Li x 

by 

wfo) = iiiyw + ff(»?(o) [/(^ ,x -') - /w] + sw* - o)[/(*f- M ) - /(»?)]. 

The operator Li x corresponds to the environment generator of a zero-range dynam- 
ics on which particles can jump between adjacent sites on each box, and between 
endpoints I and x — I. Since x < eN, we see, by adding and subtracting at most 
eN terms, that 

(f(-L l;X f)) p <(l + eN)(f(-L N f)) p . 
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Then, it is enough to prove that 



sup{({if(V(0))-#(CHO))}/)^ r - ±-(J](-LuJ]))^} 



vanishes as e J. and ( | oo. In this formula, the state space is N ' , where 
A ; * = {—/,..., 3/ + 1}, the configurations of this space are denote by the pair 
P = (^jO) where -q belongs to and ( belongs to Nq 1+1 '"' ,3!+1 ^ expectation 
is taken with respect to the measure v p l , the projection of v p on A*, Lu is the 
generator of the environment restricted to the set A^: 

(Luf)(P) = £ p(y-x)g(p(x))[f(P°«)-f(j3)] 



Epw g(m) u(P°n - m\ ; 



and the supremum is carried over all densities / with respect to v p l . 



Following the proof of the one-block Lemma 16 . 31 we need only to prove that 
supsup{ J {H(7 ] l (0))-H(C l (0))}fd^, k -^- J v7(-Li,iV?) dv Af , k } 

vanishes with limits on e and I where the supremum is on densities / with respect to 
the canonical measure fAt,k(') — v p l ("I^A* ,&), defined similarly as in (|6.1[) . where 

£ Ar , fc = {J3 e <* : (3(0) > 1, Ey eA? P(V) = k }- 

Let W0) = H(r] l (0)) - H{( 1 {0)). By the Rayleigh expansion 6, Theorem 
A3. 1.1], spectral gap estimate Lemma 15721 applied to Lij (which can be thought of 
as the environment generator on a block of length 2|A;|) and boundedness of H, 
for large N and small e, 

Wifdu At , k -if-J \fii-Lu\ff) duA t ,k 

< J W + ^Jg^ / W l {{-L l , l )-'W l }d VKhk 

< I W,<ft/ Af ,* + C(JT)J a 7e 

for some finite constant C(H) depending on H . The last term vanishes as e { 0, 
while the first term vanishes uniformly in k as I | oo by Lemma 16.61 below. □ 



Lemma 6.6. For a bounded, Lipschitz function H : R + — > R, we have that 



lim sup sup 

l^oo k>Q 



E "AT,k 



H( V l (0))-H(( l (0)) 



= 



Proof. Fix e > 0. Using that H is Lipschitz, we have that \H(r) l (0)) - H(( l (0))\ < 
C(H){r) l (0) + £ ( (0)}, and so the expectation appearing in the statement of the 
lemma is less than or equal to C(H)E UA , ^[^(O) + C'(0)]. A computation, similar 
to the one presented in the proof of Lemma 16. 4i shows that 
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for all y 7^ 0. In this formula, fJ,A*,k stands for the canonical measure defined by 

MA?,fc(*) = A*p' ('I S^gA* fii x ) = where /ip 1 is the product measure restricted 
to the set A ; *. In particular, the expectation appearing in the statement of the 
lemma is less than or equal to C{H)a 2 {\ + fc}/2|A;|. This concludes the proof of 
the lemma in the case /c/2|A/| < e. 

Assume now that fc/2|A;| > e. By definition of the canonical measure VA.*,k and 

the grand-canonical measure v p l , the expectation appearing in the statement of 
the lemma is equal to 



1 



E, 



(3(0){H(v l (0))-H(C l (0))} 



(6.6) 



Since the measure is space homogeneous, the denominator is equal to pi jk = fc/2|A;|, 
while in the numerator we may replace /3(0) by rj l (0). The numerator can therefore 
be rewritten as 



{r?(0)- Pl , k }{H(r, l (0))-H(( l (0))} 



PLkEu., 



H(r> l (0))-H(( l (0)) 



The second term vanishes because the measure /M*,fc is space homogeneous, while 
the first one is absolutely bounded by C(H)E^ A , k [ |t/(0) - pi lk \ ]. By Corollary 
6.1 (C)], this expression is less than or equal to 



C'(H)E 



1^(0) -pi.jbll < C'{H)a{pi,k)l~ 1/2 



where a(p) stands for the variance of £(0) under p p . By [5J (5.2)] and since ip(p)/p is 
bounded below and above, <j(pi i k) 2 < Cpz,fc- Therefore, if we recall the denominator 
in (16.61). we obtain that 



E„ 



H( V l (0))-H(C l (0)) 



< 



C(H) ^ 
\/lpl,k ' 



which concludes the proof of the lemma since we assumed the density to be bounded 
below by e. □ 

6.3. Proof of Proposition mi Recall H(p) = E Vp [h], Hi(rj) = H(r) l (0)), and 
Hi{p) = E Pp [H{\. Then, we have that 



E 



N 



<.N 



x=l 



< E 



E 



N 



N 



{h(t},) - H( V l s (0))}d S 

rI 1 eN 

I {H(r, l s (0))-—Y,H(v l s(v))}ds 

! 1 x—l 
t 1 eN 



E 



As h is bounded, Lipschitz, we have H is bounded, Lipschitz by Lemma 16.71 
below, and so the first and second terms vanish by Lemmatta l6.3l and l6.5l For the 
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third term, we can rewrite it as 



^{jj £ <«(x/N)(H(vl(x))-S l (ri S s N (x))}ds 



where t e (-) = e _1 l{(0, e]}. In fact, as h is bounded, we can replace t e in the 
last expression by a smooth approximation. Then, for fixed e > and I > 1, 
treating Hi(rf) — H(r] l (0)) as a local function, which is also bounded, Lipschitz as 
H is bounded, Lipschitz, the third term vanishes using Proposition 15.11 by taking 
N t oo, and e | 0. □ 

Lemma 6.7. Let h : Qn — > K 6e a local, Lipschitz function. Then, H : R + — > R 
given by H(p) — E v [h] is also Lipschitz. 

Proof. The proof is similar to that of Corollary II. 3. 7 [6] which shows h(p) = E^ [h] 
is Lipschitz. Following the proof of Corollary II. 3. 7 [6], it is not difficult to show 
{ v p '■ P > 0} is a stochastically increasing family, and for p\ < p2 that 

\H( P1 ) - H(p 2 )\ < C h Y,\E^[v{x)]-E u J V (x)]\ 

where Ch is the Lipschitz constant of h, and A c Z corresponds to the support of 
/i. If A does not contain the origin, the proof is the same as for Corollary II. 3. 7 [6]. 

Otherwise, it is enough to estimate the difference \E V [r}(fS)]—E Vp2 [r/(0)]|. When 
= p2 < pi , the difference equals E flp [?7(0) (»7(0) — l)]/pi < a<p(pi) <a 2 pia,sa~ 1 < 
g{k)/k, V {p)/p < a through (LG), (M), and E^g^Q))/^)} = <p{p)E Mp [f{ri + 1> )] 
where £>o is the configuration with exactly one particle at the origin. When < 
Pi < pi, the difference equals 

|pr% P1 lv(0) 2 } - ft % w lv(0) 2 }\ < W 2 (pi)/Pi - o- 2 (p 2 )/ Pl \ + \ Pl - p 2 \ 

where <r 2 (p) = E^ [(77(0) — p) 2 ]. The Lipschitz estimate now follows by calculating 
a uniform bound on the derivative 

q ^(P) = m 3{p) _ ° 2 {P) 

9 P Po 2 {p) P 2 

where m^{p) = Efj, p [(rj(0) — p) 3 ]. For p large, under assumptions (LG), (M), this is 
on order 0(p^ 1 ^ 2 ) from Lemma 5.2 [S] and bound a -1 < (p(p)/p < a; on the other 
hand, as p J. 0, the derivative is also bounded. □ 
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